Glueball properties in anisotropic SU(3) lattice QCD with 
improved action 
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We study the glueballs properties at finite temperature using SU(3) lattice QCD 
at the quenched level with the anisotropic lattice. We use the tree-level Symanzik 
0(a 2 ) improved action. We present our preliminary results which shows the slight 
reduction of the scalar glueball mass near T c . 



1 Introduction 

The SU (3) lattice QCD Monte Carlo simulations suggest the existence of the fi- 
nite temperature deconfinement transition at the critical temperature T c ~ 260 
MeV at the quenched levelEI. At the quenched level, without involving dynami- 
cal quark excitations, the lightest physical excitation is the scalar glueball with 
the mass 1.5 GeV < tog ~ 1-7 GeVElu, which can be expected to dominate the 
thermodynamic properties below T c . The closed packing model, a phenomeno- 
logical model for the QCD phase transition, provides us with an intuitive pic- 
ture that the deconfinement transition takes place when an "hadrons" excited 
from the thermal system begin overlap significantly with one another. How- 
ever, since the contributions from the glueballs are suppressed by a rather 
strong statistical factor as e~ mG / T , we may wonder if the total number of the 
glueballs to be excited may be insufficient in the quenched QCD. To keep the 
consistency with the closed packing model, the resolution may be found in 
either a tremendous expansion of the glueball size or a drastic reduction of the 
mass near the critical temperature (See Fig.[l]). Indeed, the lattice result of 
the correlator of Polyakov loopis. .supports the change in the shape of the qq 
potential at finite temperature UU. As a consequence, the phenomenological 



a The reduction of en glueball mass is also anticipated in the monopole condensation picture 
of the confinement □. 
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Figure 1: The glueball radius R(T C ) so as to reproduce T c ~ 260 MeV in the closed packing 
model is plotted against the effective thermal glueball mass mg(T c ). 

potential model suggests also the change in the wave function of the glueball 
particularly its mass and size. It is thus important to study the properties of 
glueballs at finite temperature, although no such studies have been done so 
far on the lattice. Even for typical hadrons, the most studies of the mass shift 
on the lattice at finite temperature are unfortunately limited to the screen- 
ing mass , i.e., the correlation along the spatial direction 0, whereas what is 
really interesting is the direct informations of the mass from the correlations 
along the temporal direction. This is because the lattice points available in 
the temporal direction decrease as the temperature becomes higher, and the 
extraction of the mass spectrum becomes more difficult. Recently, the use of 
the anisotropic lattice, whepi the temporal spacing is finer than the spatial 
one, has been established BE9, which enables us to overcome the difficulties 
mentioned above. In this paper, after analytical consideration with the closed 
packing model, we will outline how to study the glueball properties at finite 
temperature by means of the anisotropic lattice and the smearing method. 
We will present preliminary results in SU (3) lattice QCD with the tree-level 
Symanzik 0(a 2 ) improved action. 

2 The closed packing model for glueballs in the quenched QCD 

In this section, we analytically investigate the deconfinement phase transition 
at the quenched level in terms of the low-lying scalar glueball properties. Con- 
sider the bag model picture of the scalar glueball. In this picture, a sphere of 
radius R distinguishes the outside from the inside, where the gluonic degrees 
of freedom appears. As long as the temperature T is low, the number N(T) 
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of thermally excited bags is small. However, as the temperature increases, the 
number N(T) of bags increases. Finally, these bags begin to overlap each other 
and cover the whole space, and thus the deconfinement phase transition takes 
place in the closed packing model. The ratio of the volume inside the gluonic 
bags to the space volume V is given as 

_ fR(Tf x N(T) 4vr f d*k 1 m 

MT) - -R(T) J — (1) 

Now we simplify the dispersion relation to the quasi-free relation E(k,T) = 



k 2 + ttiq (T) with rriQ (T) begin the effective glueball mass at temperature 
T. Since the deconfinement phase transition is expected to take place when 
rv(T) ~ 1, the glueball radius R{T C ) near the deconfinement phase transition 
is estimated as 

-1/3 

(2) 
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(2tt) 3 eV /fc2+ m 2,(T c )/T c _ 1 J 

In Fig{t], the glueball radius R(T C ) so as to reproduce the deconfinement phase 
transition at T c = 260 MeV is plotted as the function of the glueball mass 
tog(T'c) in the closed packing model. Thus, if tog persists to be 1.7 GeV near 
T c , the glueball should have the abnormally huge radius 3.8 fm. Instead, if 
R remains to be a typical hadron size about 1 fm, the glueball mass should 
reduce drastically as tog(T' c )<0.5 GeV. 

3 The smearing method in SU(S) lattice QCD 

Using SU(3) lattice QCD, we study the glueball correlator 

G(4) = <6(t)d(0)), 6(4) = 0(t) - (O), 0(t)^Y,0{t,x). (3) 

X 

The subtraction of the vacuum expectation value is necessary for the scalar 
glueball with J PC = ++ . The summation over x physically means the mo- 
m^i^ |u|ja^ projcction as p = 0. In the case of the scalar glueball, the summand 

0(4, x) = Oia (4, x) + 23 (4, x) + 31 (4, x), (4) 

where Oij(t, x) is the real part of the trace of the plaquette operator which lies 
within the i-j plain. With the spectral representation, G(t) is expressed as as 



G(t)/G(0) = ^ Cne - B " 4 , c„= (n|O|0) /G(0), G(Q) = £ |(n|O|0) 
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where E n denotes the energy of the n-th excited state |n). (In this paper, |0) 
denotes the vacuum, and |1) denotes the ground-state ++ glueball.) Note that 
c n is a positive number with ^c„ = 1, and thus G(t) / G(0) can be expressed 
as a weighted average of the exponentials e~ Ent with the weight c„. On a 
finite lattice with the spacing a, the simple plaquette operator Oij{t,x) has 
a small overlap with the glueball ground state \G) = |1), and the extracted 
mass from G(t) looks heavier owing to the excited-state contamination. This 
small overlap problem originates from that the plaquette operator has a smaller 
"size" of 0(a) than the physical size of the glueball. This problem becomes 
severer near the continuum limit, i.e., a ~ 0. We thus have to enhance the 
overlap by improving the glueball operator so as to have approximately the 
same size as the physical size of the glueball. One jaLtiLe-aystematic ways 
to achieve this is known to be the smearing method EfHld'EJ'GJ. The smearing 
method is expressed as the iterative replacement of the original spatial link 
variables Ui(s) by the associated fat link variables, Ui(s) € SU(3) C . The fat 
link variable Ui{s) is defined to be the SU(3) C element which maximizes 



ReTr { U'(s) 



*Ui(s) UMUiis ± j)U}{s ± i) 

m ± 



(6) 



given 


as 




: 






o 


/ 




a. x 


i\ 


: 






o 





where U-^(s) = UHs — fi), and a is a real parameter. A schematic illustration 



(7) 



Note that the summation is only over the spatial direction to avoid the non- 
locality in the temporal direction. Note also that Ui(s) holds the same gauge 
transformation properties with Ui(s). We will refer to the fat link defined in 
Eq.(||) as the first fat link. The n-th fat link is defined iteratively, i.e., it is 
obtained by repeating the procedure in Eq. (J6j) with (n — l)-th fat links. For 
the physically extended glueball operator, we use the n-the smeared operator, 
the plaquette operator which is constructed with the n-th fat link variables. 
To attribute a size to the n-th smeared operator, we use the linearization and 
the continuum approximation, and obtain the diffusion equationB as 



0_ 

dn 



4>i(n; x) = DA(j)i(n] x), D 



: spatial lattice spacing. (8) 
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Here (j)i(n;x) denotes the distribution function of the gluon field Ai(x,t) after 
the n smearing operations. The n-th smeared plaquette located at the origin 
x = physically means the Gaussian extended glueball operator with the 
distribution as 



(j>i(n;x) 



1 



exp 



(47rDn) 3 / 2 

Hence the size of the operator is estimated as 
\J (x 2 ) = V 6Dn = a. 



4Dn 



Gn 



(9) 



(10) 



So far, the smearing method is introduced to carry out the accurate mass 
measurement by maximizing the overlap. However, it can be also used to 
give a rough estimate of the physical glueball size. In fact, once we obtain 
the maximum overlap with some n and a, the physical size of the glueball is 
estimated with Eq. (|Io|) . 



4 The SU(3) lattice QCD result for glueball mass and size at T ^ 

In this section, we present our preliminary results. To generate the gauge field 
configurations, we use the tree-level Symanzik 0(a 2 ) improved actional 



s G =p- 



7G 



s,i<j<3 



[c u (l - Pij(s)) + c 12 (2 - %(s) - %(«))] (11) 



+/? 7G J2 [ c ii(! - P «(s)) + ci 2 (2 - Ri 4 (s) - R M (s)) 



,i<3 



where en = 5/3, ci 2 = —1/12, j3 = 4.56 and 7g = 3.45. P^u{s) is the 
plaquette and the rectangular loop i?^^(s) is defined as i?^^(s) = ^ReTr ( 

Up(s)Up(s + pb)U v {s + 2fi)Ul{s + v + fi)Ul{s + v)Ul{s) ^. The numerical 

calculations are performed on the lattice of the sizes 16 2 x 24 x N t , N t = 28 
(T = 0.87T C ) and N t = 96 (T ~ 0)eI19. These parameters generate the spatial 
lattice spacing as a s = 0.12 fm and the temporal lattice spacing as at = 0.03 
fm. We have used 500 gauge field configurations for Nt = 28 and 150 for 
N t — 96. Throughout this section, we take the smearing parameter as a = 2.3. 

In Fig.|, the the overlap C = (G\O\0) /G(0) and x7^df are plot- 
ted against the iteration number of the smearing iV SBlear (= n). The over- 
lap becomes maximum with reasonable \ 2 as X 2 /-^df^1 in the region 10 < 
-^smear < 24. In this region, as is shown in FigM the glueball mass, which is 
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Figure 2: The overlap C is plotted against TVsmcar at T = (left figure). x 2 / N DF is plotted 
against iVsmcar at T = (right figure). 




Figure 3: The glueball mass tuq is plotted against JV smeaI (left figure). The normalized 
glueball correlator G(t)/G(0) at T = is plotted at N BrncSlT = 23 (right figure). The curve 
denotes the best single-exponential fit from the region between the vertical dotted lines. 



extracted with the single-exponential fit, takes its minimum as juq = 1.49 GeV 
at iV smcar = 23. This indicates the achievement of the ground-state glueball. 
We thus obtain mc = 1.49 GeV at T = Q. The glueball size R is estimated 
with Eq.(|To|) as 0.37 fm <i?<0.57 fm from the maximal overlap region as 
10 < -/V smoar < 24. We show also the glueball correlator G(t) at T — for 
a smearing number iV smcar = 23 in Fig.^J. Note that, if we used the isotropic 
lattice, we could only use the points indicated by the diamond symbols, which 
would be inefficient to get the mass. 

At finite temperature, we use the hyperbolic cosine fitting as G(i)/G(0) = 
C(e~ mctat + e - m G(N t -t)a t ^ where j\r t denotes the number of the temporal 
lattice size, and at the temporal lattice spacing. In Fig.[|, the overlap C and 
X 2 /-/Vdf for T = 0.87T C are plotted against A^mcar- In the whole region, we 
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Figure 4: The overlap C is plotted against iVsmear at T = 0.87T C (left figure). \ 2 /Nuf is 
plotted against iV smear at T = 0.87T C (right figure). 
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Fi gure 5: The fitted glueball mass at T — 0.87Xc is plotted against iV sm ear 

(left figure). The 

normalized glueball correator G(t) /G(0) at T = 0.87T C is plotted for Af smcar = 23 (right 
figure). The curve denotes the best hyperbolic cosine fit with the interval indicated by the 
vertical dotted lines. 



have reasonable x 2 as X /-^Vdf~1- The overlap becomes maximum in the 
region 10 < N smcai < 23. In this region as is shown in Fig|| the glueball 
mass becomes minimum as wq = 1-39 GeV in the interval 18 < N smeax < 23. 
We thus obtain m G = 1.39 GeV at T = 0.87T C . The glueball size R is also 
estimated with Eq.@ as 0.37 fm <i?<0.56 fin at T = 0.87T C from the 
maximal overlap region as 10 < iV smear < 23. We show also the glueball 
correlator G{t) at T = 0.87T C for N ameal = 23 in Fig.|. 

5 Summary and discussions 

We have applied the closed packing model to the quenched QCD to derive 
the relation between the mass and the size of the glueball so as to reproduce 
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the deconfinement phase transition at T c = 260 MeV. We have then studied 
the mass and the size of the scalar glucball at finite temperature by using 
SU(3) lattice QCD at the quenched level with the anisotropic lattice and the 
smearing method. To generate the gauge field configurations, we have used the 
tree- level Symanzik 0(a 2 ) improved action with the lattice size 16 2 x 24 xJV f , 
N t = 96(T = 0) and N t = 28(T = 0.87T C ). As for the glueball size R(T), 
no drastic change has been observed between 0.37 fm ^R(T = 0)^0.57 fm 
and 0.37 fm < R(T = 0.87T C ) < 0.56 fm. We have obtained the glucball mass 
m G = 1-49 GcV at T = and to g = 1.39 GeV at T = 0.87T C . This result 
may suggest the slight reduction of the glueball mass at finite temperature, 
which would provide a support for our conjecture discussed in Sect. 1 and 
Sect. 2. However, this reduction does not seem so drastic as is suggested by 
the closed packing model for the quenched QCD. At any rate, to draw a definite 
conclusion, much more statistics are necessary. 
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